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Abstract. We present upper limits on the amplitude of gravitational waves from 28 
isolated pulsars using data from the second science run of LIGO. The results are also 
expressed as a constraint on the pulsars' equatorial ellipticities. We discuss a new way 
of presenting such ellipticity upper limits that takes account of the uncertainties of the 
pulsar moment of inertia. We also extend our previous method to search for known 
pulsars in binary systems, of which there are about 80 in the sensitive frequency range 
of LIGO and GEO 600. 
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1. Introduction 

Over the last few years the LIGO and GEO 600 gravitational wave detectors have had 
three (SI, S2 and S3) and two (SI and S3) science runs respectively. Data from SI have 
been used to set upper limits on gravitational wave emission from a variety of possible 
sources PI El El E] • Here we concentrate on a search for gravitational waves from known 
pulsars which has been extended both in terms of analysis method and number of sources 
from that presented in pp. 

SI (23 August - 9 September 2002) was the first major data taking period and 
a search for gravitational waves from only one pulsar (the fastest millisecond pulsar 
J1939+2134) was performed [TJ. This search used both a time domain and a frequency 
domain analysis method. For S2 (14 February - 14 April 2003) the number of pulsars 
being searched for increased to 28 [S], with the criterion that the pulsars were isolated 
and had spin frequencies within the sensitive region of the instruments, i.e. / > 25 Hz 
or gravitational wave frequencies of few > 50 Hz, assuming emission from a triaxial 
neutron star [6 . The time domain analysis technique was considered the optimal search 
for S2 where the source position and frequency parameters were known. The analysis 
and results from S2, including possible astrophysical interpretations, will be summarised 
in this paper. 

Analysis of the S3 run (31 October 2003 - 9 January 2004) is underway and has 
the aim of extending the number of pulsars in the search to all those with / > 25 Hz 
(~ 110) including those within binary systems. The additions to the analysis needed 
to incorporate the binary systems will be discussed in more detail in this paper. We 
will also discuss how the results on neutron star ellipticity could be interpreted as an 
exclusion region on a moment of inertia - ellipticity plane. 

2. S2 summary 

2.1. Analysis 

S2 data were available for analysis from all three LIGO interferometers. Frequency and 
positional information for all the isolated pulsars with / > 25 Hz can be found on the 
Australia Telescope National Facility (ATNF) ^ online catalogue. However for some 
pulsars this information lacked the accuracy required to be sure that the search was 
coherent over the period of the science run, and additional timing over the epoch of S2 
was requested and obtained for 18| of these objects |Hj. In total this gave us 28 pulsars 
with reliable information upon which to reconstruct their phase evolution. Of these 28 
pulsars, 14 are in globular clusters. The list also includes the Crab pulsar (J0534+2200) 
and the fastest millisecond pulsar J1939+2134 targeted in SI. 

The analysis method is described in detail in [5]. In summary we perform a time 
domain heterodyne of the data with the known phase evolution of the pulsar signal. 

t The Crab pulsar (J0534+2200) data was taken from the onhne Jodrell Bank Monthly cphemeris 
,http : //www . jb . man . ac . uk/research/pulsar/crab . html 
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no. of pulsars 


1x10-24 < ho < 5x10-24 


20 


5x10-24 <ho< 1x10-23 


4 


ho> Ix 10-23 


4 



Table 1. The 95% upper limits on ho for the 28 pulsar searched for with the S2 run. 



ellipticity e 


no. of pulsars 


1x10-6 < e < 1x10-5 


4 


1x10-5 < e < 1x10-4 


16 


e > 1x10-4 


8 



Table 2. The upper limits on e for the 28 pulsar searched for with the S2 run. 



After heterodyning, any signal in the data would vary only with the beam pattern [Hj of 
the detector and its form would depend on the four unknown parameters of gravitational 
wave amplitude (/iq) and polarisation angle (tp), pulsar orientation angle (t) and the 
signal's initial phase (0o)- The heterodyned data is then filtered and re-binned from 
the detector sample rate of 16384 Hz to 1/60 Hz. We then determine the probability 
distribution functions of the parameter values using a Bayesian inference technique, as 
shown in [Hj. 



2.2. Hardware injections 

During S2 we injected two artificial pulsar signals into the three LIGO interferometers 
for 12 hours by modulating one mirror of each via the actuation control signal. This was 
done to perform a validation of the search pipeline from as far up the chain as possible. 
The injections increased our confidence in the phase calibration of the detectors, and 
allowed for a joint coherent analysis of data from all the detectors. 



2.3. Results and astrophysical interpretation 

The 95% upper limits on ho are summarised in Table. ^ We can recast this upper limit 
on ho as an upper hmit on the pulsar's ellipticity, e, calculated via 

ho r 1 Hz2l038kgm2 

where r is the distance to the pulsar in kpc, / is its spin frequency, and I^z is its 
principal moment of inertia jHj . The ellipticity limits are summarised in Table |21 using 
the canonical value of 10^^ kg m2 for the moment of inertia and assuming no error on the 
pulsar's distance. The lowest value of strain upper limits (1.7x10-24) is associated with 
pulsar J1910-5959D and of ellipticity (4.5 x 10-^) with pulsar J2124-3358. The upper 
limit for the Crab pulsar, at ho = 4.1 x 10-2^, is only a factor of ~ 30 greater than the 
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canonical limit inferred from simple spin-down arguments§. For all other pulsars the 
upper limits are a few orders of magnitude higher than their canonical limits. 

Although the inferred ellipticities are well above those allowed by the majority 
of conventional neutron star equations of state they are approaching the region of 
astrophysical interest. Indeed, for the lowest pulsar ellipticities (~ 4.5 x 10~^) we are 
starting to reach the range permitted by at least one exotic theory of neutron star 
structure JT] . It should be stressed that these results are the first direct upper limits 
on gravitational wave emission for 26 of the 28 pulsars, with previous upper limits for 
pulsar J1939+2134 and the Crab pulsar given by [T] and [TU] respectively. For the pulsars 
within globular clusters the spin-down can be masked by local Doppler shifts caused by 
the cluster dynamics. This makes placing a spin-down based upper limit difficult unless 
the pulsar motions within the clusters can be found independently. These gravitational 
wave interferometer upper limits are inherently independent of the cluster dynamics. 

3. S3 analysis 

We expanded the search in the S3 analysis to include all pulsars with / > 25 Hz. These 
include ~ 70 pulsars within binary (or more complex) systems so it is important to 
account for the extra Doppler and relativistic time delays due to the pulsar's orbital 
motion. An advantage over S2 is the availability of GEO 600 data in addition to LIGO 
data. 

3.1. Binary pulsar signal 

In the case of an isolated pulsar the signal received at the detector needs to be corrected 
to the solar system barycentre (SSB) by calculating the Doppler delays and other 
relativistic effects. This is possible as the pulsar's position is known and we have very 
good solar system ephemerides. Of course it is simply the reverse procedure to that 
used to determine the positional and frequency parameters of the pulsar from radio 
data, so the reconstruction is generally of good quality. The motion of a pulsar in a 
binary system adds a number of Doppler and relativistic time delays: 

AT,i„ = Ai, + Ae + As + Aa, (2) 

where Ar is the Roemer delay (light travel time), Ae is the Einstein delay due to 
special relativistic effects, Ag is the general relativistic Shapiro delay and A^ is the 
abberation delay caused by the pulsar's rotation. These delay changes can be far more 
pronounced than those from the Earth's motion, with up to a 0.03 Hz frequency shift. 
The time delays are parameterised by properties of the binary system including its 
period, eccentricity, angular velocity, time of periastron, projected semi-major axis and 
several relativistic parameters depending on the nature of the system. These parameters 
are found by fitting the radio observations (using the standard TEMPO data reduction 

§ This limit is imposed by energy conservation and assumes that the rotational kinetic energy is only 
lost to gravitational waves 
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package [12 ) to various binary models (e.g. The model to which the binary system 

is fitted will depend on how relativistic it is or which parameters you wish to extract. 
The 70 binary system pulsars fall mainly into two models: 32 in the low eccentricity 
(ELLl model and 33 in the Blandford-Teukolsky (BT model JHI)- Four fall into 
the highly relativistic D amour- Deruelle (DD ^3]) naodel and a further one into the 
Blanford-Teukolsky-2-Planet (BT2P) model, although this can be adequately fit using 
the simpler BT model. With this additional binary information provided the full set of 
pulsars can be included in the search. So far new timing information for the majority 
of these pulsars has been provided [H] and the analysis is underway. 



3.2. Hardware injections 

We injected ten artificial isolated pulsar signals, with a wide range of signal parameters 
into the LIGO interferometers. One other signal was injected into both GEO 600 and 
LIGO. The signal strengths ranged from marginally detectable to very strong. These 
signals have been successfully extracted using the search pipeline jH]. For the strongest 
signals, agreement between the injected and recovered parameter values has been to 
within a few percent. 



4. Moment of inertia - ellipticity plane 

So far we have used the canonical value for the neutron star moment of inertia 
(lO^^kgm^) when calculating the pulsar ellipticities. This is the moment of inertia of a 
1.4 Mq sphere of uniform density and radius 10 km. The true value can vary by factors 
of a few under different models for the neutron star equation of state. Equation [H shows 
that Hq can be used to set an upper limit on the neutron star's quadrupole moment, 
Izz^, that is independent of the actual value of Izz- This value can then be used to form 
an exclusion region in the Izz — e plane (see Fig. Q). This plane can provide exclusion 
regions on both the moment of inertia and ellipticity or can be used to read off an upper 
limit on ellipticity for a preferred value of Izz such as the canonical value. 

Figure |21 shows such an exclusion region for the Crab pulsar upper limit from S2. 
An exclusion region can also be set using the spin-down based upper limit argument, at 

Z2G{2nffe- ^' 
The point at which the experimental results and the spin-down result cross shows the 
the point at which the experimental results beats that of spin-down. For the S2 Crab 
pulsar result it can be seen that the experimental result only beats the spin-down result 
for unfeasably large values of Izz (> 8xl0^°kgm^). The upper limit on the ellipticity 
for the experimental result and assuming Izz = lO'^^kgm^ is ~ 2x 10~^. 
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Figure 1. This figure shows a rough representation of the regions in the moment 
of inertia I^z - ellipticity e plane that can be excluded via various methods. The 
electromagnetic emission of a pulsar can set a lower limit on the moment of inertia 
by equating the EM emission with the rotational energy loss of the pulsar The 
various equations of state for neutron stars can constrain the mass/radius relation and 
therefore moment of inertia. Equations of state will also put limits on the maximum 
allowable ellipticity of the neutron star. A limit can be set from upper limits on 
gravitational wave emission. 




Figure 2. The areas above the full lines are the excluded regions on the moment of 
inertia - ellipticity plane for the Crab pulsar as obtained from the S2 ho upper limit 
and the spin-down arguments. These assume a distance to the Crab pulsar of 2 kpc. 
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5. Future work 

The analysis of S3 data is presently underway targeting most known pulsars within 
the band. In addition we will be incorporating other sources of error into the current 
Bayesian framework, including systematic calibration errors in the value of ho and pulsar 
distance errors into the value of e or the I^z — e plane exclusion region. The S3 results 
promise to give a reasonable improvement over S2 for some pulsars. In particular, the 
Crab pulsar result should only be a factor of a few above the spin-down limit. 

The next science run of LIGO and GEO 600 (S4) is due to be underway in early 
2005 providing more data to analyse. The pipeline used to obtain the results in this 
paper is at a mature stage and will be applied directly to the S4 data. 
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